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Abstract 

Research on the interplay between the dynamics on the network and the dynamics of the network has attracted much 
attention in recent years. In this work, we propose an information-driven adaptive model, where disease and disease 
information can evolve simultaneously. For the information-driven adaptive process, susceptible (infected) individuals 
who have abilities to recognize the disease would break the links of their infected (susceptible) neighbors to prevent the 
epidemic from further spreading. Simulation results and numerical analyses based on the pairwise approach indicate 
that the information-driven adaptive process can not only slow down the speed of epidemic spreading, but can also 
diminish the epidemic prevalence at the final state significantly. In addition, the disease spreading and information 
diffusion pattern on the lattice give a visual representation about how the disease is trapped into an isolated field 
with the information-driven adaptive process. Furthermore, we perform the local bifurcation analysis on four types of 
dynamical regions, including healthy, oscillatory, bistable and endemic, to understand the evolution of the observed 
dynamical behaviors. This work may shed some lights on understanding how information affects human activities on 
responding to epidemic spreading. 
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1. Introduction 


The spreading dynamic is one of the core issues in network science |[ll0]. where most of the related research 
processes focus on epidemic spreading and information diffusion in recent years. Much of the work to date focuses on 
the analysis of these two processes independently, such as the spread of single contagion 0 or concurrent diseases 
|2l, and the diffusion of various kinds of information (e.g., news |Q], rumor 101, innovation |0].). However, the epidemic 
spreading process is closely coupled with the corresponding disease information diffusion (or saying awareness) in the 
real world. For example, the Google Flu Trend (GFT) can detect regional influenza outbreaks with the online search 
for the disease-related information fid. During the severe acute respiratory syndrome (SARS) outbreak in China 
in 2003, overwhelming number of disease reports have helped people win the ’’people’s war against SARS” Anally 
idd. It is obvious that disease information diffusion plays an important role in control of the epidemic outbreak, 
but it is not easy to quantitatively measure the strength of its impact lll3n . 

Nowadays, some modelin g approa ches on the interaction between epidemic spreading and disease information 
diffusion have been proposed |3 14, 15]. The fundamental assumption is that, when a disease breaks out, people 
may get the disease information from their friends or media before the advent of the epidemic and take some pre¬ 
ventive measures to keep away from being infected 114, dd. By depicting preventive measures as the reduction 
of transmitting probability |18, 19] or particular states of individuals (immune or vaccination) Idl . previous mod¬ 
els showed that the disease information diffusion indeed inhibits the epidemic spreading significantly (reduce the 
epidemic prevalence and enhance the epidemic threshold) Cl EH]. Therefore, the emergence of mutual feedback 
between information diffusion and epidemic spreading lll3n exhibits the intricate interplay between those two types 


Email addresses: liuchuang@hznu . edu . cn (Chuang Liu), zhangzike@gmail. com (Zi-Ke Zhang) 


Preprint submitted to Elsevier 


May 20, 2015 















of spreading dynamics. It should be noted that such complex interacted spreading dynamics are all performed on 
fixed networks (or population) in the aforementioned studies. However, individuals would sometimes cut off the con¬ 
nections with the infected ones when they are aware of the disease, leading to the continuous evolution of network 
structure. Consequently, how to characterize the mutual spreading process on the dynamical networks is a crucial 
issue we want to address in this work. 

Generally, the network dynamic researches could be classified into two lines; (i) one is treated as the dynamics 
of the network, which focuses on the time evolution of network topological properties 11221123112411 : (ii) the other is 
considered as the dynamics on the network, which concerns the dynamical process of the nodes (or interactions) on 
networks, such as the epidemic spreadin g an d information diffusion process |25, 2^, the evolutionary game 1271 

1281 1 becomes the most important model to describe the mutual interaction 


and so forth. Currently, adaptive process 
between temporal network topology and node dynamics. In the adaptive process of epidemic spreading dynamics, the 
susceptible individuals are allowed to protect themselves by rewiring their links from the infected neighbors to some 
other susceptible ones ^ Many researches indicate that segregating infected (or susceptible) individuals 
with the adaptive process is an efficient strategy to reduce the fraction of susceptible-infected {S /) interactions, as 
well as hindering the outbreak of the whole epidemic spreading process ^2, 33, 34]. In addition, abundant temporal 
behaviors are presented for illustrating the spreading dynamics on the adaptive network, such as the coexistence of 
multiple stable equilibrium and the appearance of an oscillatory region, which are absent in the spreading dynamics 
on static networks ll2^l35n. Besides the edge rewiring strategy, the link cutting or temporarily deactivating is also a 
commonly used rule Obi 13711 . 

In this work, we aim to understand the complicated interplay between the epidemic spreading and disease infor¬ 
mation diffusion on dynamical networks. We propose an information-driven adaptive process to model the dynamical 
networks, where those who have been informed of the emergence of disease can break their neighbouring connections 
to avoid further epidemic spreading. And epidemic spreading and disease information diffusion are described by the 
SI and SIS model respectively, where the disease information generation of the infected individuals is considered to 
form a mutual feedback loop between these two types of spreading dynamics 11811 . Therefore, the effect of information 
diffusion on epidemic spreading could be interpreted by two aspects: (i) reduce the epidemic spreading probability 
with protective measures; and (ii) cut off SI links with the information-driven adaptive process. Both numerical anal¬ 
yses based on the pairwise approach and simulation results indicate that the information diffusion and the adaptive 
process can inhibit the epidemic outbreak significantly. In addition, we present a full local bifurcation diagram to 
show the complex dynamical behaviors in the proposed model. 


2. Results 

Model description. In Fig. [1] we illustrate the proposed model in detail. The vertical transformation describes the 
diffusion of disease information by an SIS model, where individuals can be at one of the two states: (i) -I-: indicates 
that the individuals have known the existence of the disease, denoted as the informed ones; (ii) -: indicates that 
the individuals have not known the existence of the disease. At each time step, the informed nodes will transmit 
the information to their unknown (-) neighbours with probability a, and each informed individual can forget the 
information of the disease with a probability A. Besides, the one who has been infected by the disease will become to 
know the information of the disease with a corresponding rate o) IBIIl. 

In the horizontal transformation of Fig. [T] the epidemic spreading is described by an SI model. Each node is at 
one of two states, susceptible (S) or infected (I). The disease can be transmitted through the SI links, where the S-state 
individuals could be infected with the probabilities fi, crijA, crsji and ctsiP respectively through S-I-, S^I+, S +1- and 
S +1+ links, where cri,crs and cr^/ are the impact factors of the information on epidemic spreading. Generally, when 
people know the occurrence of the disease (informed individuals), they would like to take some measures to protect 
themselves, leading to the reduction in infectivity (0 < (Ts,o-j < 1). In particular, the influence coefficient of the 
epidemic spreading probability through S +1+ links could be calculated as cr^/ = o-gcri, with the assumption of the 
independent effect of the infected probability. 

Additionally, we consider an information-driven adaptive process which the informed individuals would reduce 
physical contacts to protect themselves or their friends. That is to say, the informed susceptible individuals {S+) will 
keep away from their infected neighbors to protect themselves from being infected, and informed infected individuals 
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(/+) will also avoid contacting their susceptible neighbors to prevent the epidemic from further spreading. Conse¬ 
quently, the edge-breaking rule of adaptation is adopted Thus, at each time step, the S+ (/+) state individuals 
will break the links connected to their I (S )-state neighbors with probability rs (rj) respectively. Specially, the break¬ 
ing probability of the S +1+ pairs could be interpreted as 1 - (1 - )(1 - r/) with the independent assumption. It should 

be worth noting that the deactivation of SI links only represents the avoidance of physical contacts between the S - 
and /-state individuals. That is to say, the edge-breaking process will not affect the diffusion of disease information 
for it can be transmitted through other types of connections such as phone, internet and so forth. It is worth noting 
that the dynamic of the epidemic spreading is simply a classical SI model when the parameters are set as rs - ri - 0, 
i.e., there is no edge-breaking in this case. 

According to the model described above, the spreading process can be summarized as follows. At the beginning, 
an individual is randomly selected as the /+ node, which is considered as the seed of both the epidemic spreading and 
information diffusion, and all other individuals are set as 5 _ ones. At each time step, (i) the infected individuals would 
transmit the disease to their susceptible neighbors with the corresponding probabilities; (ii) the informed individuals 
would transmit the disease information to their un-informed neighbors; (iii) the informed individuals can forget the 
information; (iv) the informed individuals would also break the links with their relevant neighbors by considering the 
adaptive mechanism. Finally, the spreading process would be terminated when the size of the infected individuals 
becomes stable. 


Numerical mathematical analysis. Firstly, we develop theoritical analyses to depict the dynamic processes of both 
information diffusion and epidemic spreading. In particular, mean-field analysis and the pairwise analysis are adopted. 
Let X be the state variable, thus [y] denotes the expected values of individuals of different type at the population 
(e.g. [5'+] and [5'+/+] represent the expected number of informed susceptible nodes and expected number of links 
connecting an informed susceptible node to an informed infected node respectively). 

Therefore, with the classical mean-field approach, we can obtain: 


d[h] 

dt 


(k>[5+](crs;8[/_] + asCTiPiU]) + a[I.]([S^] + [/+]) + tu[/.] - A[h] 


comparatively, with the pairwise approach, we can obtain: 


d[h] 

dt 


{(TsPiS +/-] + CTso-iPiS +/+]) -H a([5 +/-] -H [/-/+]) -H u)[I-] - T[/+] 


( 1 ) 


( 2 ) 


where, the first terms of Eq.([T]) and (|2]i describe the infection of the S +-state individuals, the second terms describe the 
information acceptance of the /--state individuals, the third terms describe the information generation of the /--state 
individuals and the last terms represent the information loss of the /+-state individuals. Simultaneously, the full set of 
differential equations based on those two approaches can be illustrated in Appendix. By the way, the adaptive process 


could be described by the last terms of 


d[S+I-] ^/[5-/+] . ^/[5+/+] 

and 


dt dt dt 

be noted that the pairwise analysis is based on a well-known closure approximation given by [ABC] - 


in the pairwise approach of Eq. dUi. It should 

[AB][BC] 


[B] 


with the assumption that the degree of each individual obeys Poisson distribution 113811391] . In general, it might be very 
hard to get exact solutions of such complex differential equations, thus we give numerical solutions of the equations 
instead of the theoretical analysis in the following analysis. 


Simulation and numerical analysis without adaptive behaviour. We first consider a simple case of no adaptive 
process, i.e., there is no edge-breaking. Eig. |2]gives the simulation result of the epidemic spreading dynamics for 
various information diffusion probabilities a, with the epidemic spreading probability p - 0.3. Eor the SI process, the 
whole population would be infected when p > Q for the connected social networks, resulting in that the infected density 
equals to 1 at last for all the value of a in Eig. |2] That is to say, the diffusion of the disease information couldn’t 
reduce the size of the epidemic outbreak with the SI assumption. Whereas it is obvious to find that the epidemic 
spreading speed decreases with the increase of a, where the time cost for the whole population becomes infected 
when a = 1 is about three times longer than that of o' = 0. In this sense, the diffusion of the disease information can 
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inhibit the epidemic spreading significantly. In addition, the inset of Fig. [^indicates that the epidemic spreading can 
enhance the disease information diffusion. Actually, according to the transformation process in Fig. [T] on the one 
hand, we can realize that the epidemic spreading could be influenced by information diffusion where the epidemic 
spreading probability of the informed individuals would change; and on the other hand, the information diffusion 
could be influenced by the epidemic spread where the social disease information level (namely Info in the inset of 
Fig. |2]i would be higher if more people are infected for the information generation, denoted by the parameter u. In 
this way, a mutual feedback between disease spreading and information diffusion emerges: higher prevalence of the 
infected individuals keeps more disease information generation in the population, which in turn gives rise to more 
informed individuals, thereby weakening the spread of epidemic. 

Fig. [3 presents the numerical analysis and the simulation result of the epidemic spreading dynamics when 
rs - ri - 0, and the curves of the mean-field approach and pairwise approach are calculated according to Eq. 
Q and (IHi respectively. The epidemic spreading based on the classical mean-field approach is much quicker than 
that of the simulation result, which would be caused by the mean-field assumption on the SI model. In the mean-field 
assumption, the I- and S -state individuals are well-distributed in the system. However, in the SI process, the /-state 
individuals are all well clustered, resulting in that many /-state individuals have no chance to contact the S -state in¬ 
dividuals. In this way, the classical mean-field approach can not exactly describe the SI dynamic process. However, 
such problem is not so significant in the pairwise approach, which consider the time evolution of the links directly. 
And Fig. |3] shows that the time evolution of the spreading process on the pairwise approach finds good agreement 
with the simulation result. 

Spreading dynamics with the adaptive process. In this part, we shall present the spreading dynamics with the 
information-driven adaptive process, shown in Fig. 0] Different from the results of Fig. |2] the saturation value of the 
infected density at the final state is much smaller than 1 in Fig. 0^. That is to say, with the adaptive process based on 
the information diffusion, many individuals could avoid being infected via reducing some contacts. In addition, we 
also plot the numerical solution based on the pairwise approach in Fig. 0^. It can be seen that the pairwise solution is 
not so well for the spreading dynamic on the adaptive system. The difference might be caused by the network structure 
variation in the adaptive process, where the assumption of the pairwise approach is the Possion degree distribution. 
And this conjecture is proved in Fig. ih, where the degree distribution of the original network is approximate to 
the Possion-distribution (pink circle markers), while the distribution of the network at the final state (gray diamond 
markers) is quite different. In addition. Fig. 0^ shows that the difference becomes larger with the increase of time, 
where the degree distribution deviates more away from the Possion distribution when the process goes on. 

The information-driven adaptive process can not only slow down the speed of epidemic spreading, but can also 
diminish the epidemic prevalence at the final state significantly according to Fig. |2]and Fig. 0] For simplicity, we 
assume rs - rj - r'm the following analysis. In order to exhibit the influence of information diffusion in detail, we 
show the full phase diagram a - p for r = 0.1 in Fig. |5] The Fig. |5^ and|5]5 are the numerical solution of the pairwise 
approach and the simulation result, respectively. As stated previously, the numerical solution is not very precise, but 
it can match the overall trend of simulation result well. For a fixed epidemic spreading probability p, the size of 
the epidemic outbreak reduces with the increase of a. That is to say, the disease information diffusion can inhibit the 
epidemic spreading. Analogously, the quicker and broader of the information diffusion (larger a) is, the more efficient 
inhibition on the epidemic spreading will be. In addition, the curve of the color mutation (the dashed green curve) in 
Fig. |5]could be considered as the transition point, where the epidemic can’t spread out if a and p locate at the area on 
the left of this curve (the white range). The threshold value of the epidemic spread probability becomes larger with 
the increase of a. 

In order to intuitively demonstrate the epidemic spreading and the information diffusion process, we show the 
simulation results of those two types of spreading processes on a 100 x 100 lattice in Fig. |6]for various a. It presents 
four kinds of different levels of information spreading processes (corresponding to different a), and observe how the 
information diffusion affects the spreading of epidemic. In addition, as the adaptive edge-breaking process is merely 
executed on the epidemic spreading process, while these edges can still transmit information, the prevalence of in¬ 
formation diffusion could maintain at a high level. For each a in Fig. firstly we give the fraction of the infected 
individuals in each time step (the red curve in each subfigure). And for some particular time steps, we show the 
states of each individual with the gridding patterns, where the red dots and the gray dots represent the infected and 
informed individuals respectively (the contact networks and the un-informed susceptible individuals are not shown in 
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the figures). We can intuitively see the distribution of the infected and informed individuals and conclude that when 
the diffusion of information is slower than the epidemic, we cannot stop the epidemic from spreading (Fig. and 
Eh), however, when the information is diffusing faster, the epidemic will be trapped into an isolated area and cannot 
spread anymore (Fig. |6j; and|6}l). 

The sensitivity of the edge-hreaking prohahility on epidemic spreading dynamics. The phase diagram in Fig. 
|5] shows the impact of information diffusion rate a on the epidemic spreading dynamics. In general, the adaptive 
edge-breaking probability rs and r/ are also important parameters in affecting the epidemic spreading process. Fig. [T] 
illustrates the epidemic prevalence in the final state versus the adaptive edge-breaking rate (r) for various information 
diffusion rate a. It can be found that the epidemic prevalence diminishes with the increase of r, that is to say, the 
epidemic could be controlled if people are very sensitive with the disease information and subsequently keep away 
from the infected. It should be noted that there is no disease information diffusion when a = 0, but with considering 
the information generation, the infected individuals could stop contacting with the susceptible neighbors to impede 
the further spreading of epidemic. With the increase of a, the epidemic prevalence reduces sharply versus r and the 
continuous transition could be observed. By the way, it will change to a total isolation of infected individuals for 
r = 1, which seems to be the most effective way in controlling the contagion SEl]. 


Abundant dynamical features based on information-driven adaptive process. In order to deeply characterize the 
complex dynamical features of the proposed process, we concentrate on the distribution of the infected density in the 
final state (/*) rather than the simple average value |28, Fig. 0shows four different types of dynamical behavior 
by calculating the distribution of the final fraction of infected for various p and r . For the distribution of Fig. [^, we 
have carried out 10000 realizations of the infected density, and above 94% of the infected density is 0.0001, and the 
maximal is 0.0007, i.e., the infected density I* —> 0, thus we consider this distribution indicates a healthy state (the 
disease can’t spread out) under the parameters setting here. Similarly, as to the distribution of Fig. 01, above 90% 
of the infected density is higher than 0.8, indicates a case of endemic state (epidemic outbreaks). Whereas the case 
illustrated in Fig. 0: is very different, where the infected density /* is around either zero or a nonzero value. This 
indicates that a bistable state 11281] is located in this model, where healthy state and endemic state are both stable in 
this case. In addition, the oscillatory dynamic can also be observed in particular parameter settings (Fig. 03). 

According to the dynamical behavior illustrated in Fig. 0under different parameter sets. Bifurcation diagram of 
the density of the infected as a function of infected probability fi for different values of the edge-breaking rate r is 
given by Fig. |9}i. Without the adaptive edge-breaking mechanism (r = 0), the disease can spread out only if j3 > 0 
for the SI process. When r > 0, the dynamical behaviors become more complicated, where the discontinuous phase 
transitions, bistable, oscillatory are observed. A fast edge-breaking (large r) leads to a broad healthy and bistable state 
range (shows by the range in the arrows) in Fig. Ell. In Fig. Eb, we give a full r - p bifurcation diagram according 
to our simulation results, and we can clearly identify the areas of healthy, oscillatory, bistability and endemic state 
in this model. At last, we present the dependence of the average value of infected density over 10,000 independent 
realizations on r and p in Fig. Et, where the changing of the density is consistent with the area classification in Fig. 

Wp- 


3. Conclusion 

In order to understand the interplay between the dynamics on the network (the spread of epidemic spreading 
and disease information) and the dynamics of the network (the time varying of network links), we present two types 
of spreading dynamics with SI and SIS process respectively on an information-driven adaptive network, where the 
individuals who have known the disease information would probably cut off their links with others. Firstly, we 
illustrate the mutual feedback of epidemic spreading and information diffusion without considering the edge-breaking 
process (rg — rj - 0), where the high epidemic prevalence preserves high disease information level, which in turn 
slows down the epidemic spreading. In this case, the numerical analysis based on the pairwise approach is consistent 
with the simulation result very well. Secondly, the results are very different when the information-driven edge¬ 
breaking process is considered (r^, r/ > 0). The epidemic cannot spread out if the spreading probability is smaller 
than the threshold (shown in Fig. 0). In addition, the disease spreading and information diffusion pattern on the lattice 
give a visual representation that the disease might be trapped into an isolated field with information-driven adaptive 
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process. Therefore, the information-driven adaptive process can inhibit the epidemic spreading significantly that it 
can not only slow down the epidemic spreading speed, but also reduce the epidemic prevalence. Finally, we give the 
local bifurcation analysis on four types of dynamical phenomena, including healthy, oscillatory, bistable and endemic, 
indicating that the state changes from healthy to oscillatory, bistable, endemic state as p increases. 

In summary, we study the dependence of the epidemic spreading on the disease information diffusion and the 
information-driven adaptive process, with considering the simplest spreading model (SI) and adaptive proces s (edge - 
breaking). Recent researches show the different features between the epidemic and the information diffusion l42ll43n . 
and this difference dynamics would also impact the interplay between epidemic spreading and disease information 
diffusion significantly. Another area for future extension is to adopt other adaptation rules rather than the simple 
edge-breaking strategy, such as the temporarily deactivating, where the broken links would be active again after a 
fixed time 1371 or if the corresponding infected node becomes recovered | 


Acknowledgments 

This work was partially supported by Natural Science Foundation of China (Grant Nos. 11305043 and 11301490), 
Zhejiang Provincial Natural Science Foundation of China (Grant No. LY14A050001), Zhejiang Qianjiang Talents 
Project (QJC1302001), and the EU FP7 Grant 611272 (project GROWTHCOM). 

Appendix 

Denote [x] as the expected values of individuals of different type described in Sec. Numerical mathematical 
analysis, the epidemic spreading is depicted by the parameters fi, (JiP, o-s/3 and crsi/3, while the diffusion of disease 
information is controlled by the parameters: a. A, a>. All these parameters have been explained in Sec. Model descrip¬ 
tion. According to the model described above, the the differential equations of the mean-field approach and pairwise 
approach are given as follows. 

Mean-field approach: 

= - <k>/?[/_][5_] - {k)cr,/3[h][S.] - a([5+] + [/+])[5_] + 
at 

= - {k)crs/3[I-][S^] - {k)(rscr,/3[h][S^] + cr([5+] + [/+])[5_] - 
at 

(3) 

^ =(k)^[/-][5_] + {k)(rili[h][SA - a([5+] + [/+])[/-] - a)[IA + d[/+] 
at 

^ =(k)(rs/l[/_][5+] + {k)ascril3[I^][S^] + a([5+] + [/+])[/.] + w[/_] - A[h] 
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Pairwise approach: 

45-] 


dt 

45+] 

dt 

d[I-] 

dt 

dlh] ^ 

dt 

d[S-I-\ 

dt 


= -45-/-] - (TiPiSM] - 4[5-5+] + [5_/+]) + A[S^] 

= - cTsPlSJ.] - cTs(riP[S^U] + a([5_5+] + [5-/+]) - /l[5+] 
=^[5-/_] + (r/^[5_/+] - a([5+/_] + [/-/+]) - w[/-] + A[h] 
^(TsPiS +/-] + o-sctjPIS +/+] + ai[S +/_] + [/-/+]) + w[/-] - /![/+] 


= -45-/-] + 4[5+/-] + [5_/+]) - w[5_/-] +p 


[5-/_]([5_5_] - [5-/_]) 


[5-] 


■CTlP 


[5_/+]([5_5_]-[5-/-]) 

[5-] 


[5_/_]([5-/+] + [5_5+]) [5 _/-]([/-/+] + [5+/-]) 

■ a -——- a - 


[5-] 


[/-] 


45-/+] , n , rc , n , ,rc , n rc , n ,rc , n J5_/_][5-/+] 45-/+]^ , JS J^][S ,] 

= -CTjPIS-1+] + w[5_/_] + /i[5 +/+] - a[5_/+] - /{[5_/+] —p -j—-j- CTlP -h CTsP- 


dt 


[5-] 


[5-] 


[5+] 


45+/+][5-5+] [5-/_]([/-/+] + [5+/-]) [5_/+]([5_/+] + [5-5+]) 

+ CTsCTiP --h a -- a -- rpd _/+] 


[5+] 


[/-] 


[5-] 


d[S+I-] orrri.irrri ircri rcri rcri o [5+/+] [5+/_] 

= -crsP[S+I-] + /i[5+/+] - /l[5+/-] - Qr[5+/_] - w[5+/_] - (TsP - cTso-iP- 


dt 


[5+] 


[5+] 


45_/_][5-5+] 45-/+][5-5+] [5-/_]([5_5+] + [5_/+]) [5+/_]([/-/+] +[5+/-]) 

+ P -——- + CTjP --h a --O'-;-r5[5+/_] 


[5-] 


[5-] 


[5-] 


fi^[5+/+] 

dt 


-crs(TiP[S +/+] + a[5_/+] + a[5 +/_] + w[5 +/_] - 2A[S +/+] + osP 


[/-] 

[5+/-]([5+5+]-[5+/+]) 

[5+] 


^[5+/+]([5+5+] - [5+/+]) [5-/+]([5_/+] + [5_5+]) [5+/_]([5+/-] + [/-/+]) 

+ CTsCTiP - T^r-- -1- a - T^r-- -H a- 


[5+] 

[l-(l-r4(l-r/)][5+/+] 


[5-] 


[/-] 


= 245-/-] . 24/-/+] - 24/-/-] . 2^^ . 


dt 

</[/-/+] 

dt 


[5-] [54 

= (TiP[S-h] + asPiS^I-] + 4[/-/-] - [/-/+]) + ^([/+/+] - [/-/+]) - «[/-/+] 


[/-] 


45_/4[5-/+] J5+/42 45+/+][5+/4 [/-/-]([5+/4+ [/-/+]) 

+ P -——-+ CTiP -h (TsP + ctsctjP --h a- 


- a- 


[5-] " [54 

[/_/+ ]([5+/4 + [/-/+]) 
[/-] 


[5+] 


[5+] 


[/-] 


d[I+I+] r 1 , T rr r 1 , o rr r 1 oirr r l , T „ [5+/-] [5+/+] [5'+/+]^ 

= 2crscr/P[S +/+] + 2ci'[/-/+] + 2ct)[I-I+\ - 2A[I+I+\ + 2crsP ----H 2crs<TiP- 


dt 


[5+] 


[5+] 


+ 2Qr 


[/-/+]([5+/4 + [/-/+]) 
[/-] 


45-54 c n ,4^-'^-][5-54 ,_45-/+][5-54 ,45_54([5-5+] + [5-/+]) 

— 2/x[o _kj -I-] — 2y5 — 2cr jp ^ ^ — 2ci'“ 


dt 

45-5+] 

dt 


[54 


[5-] 


[5-] 


,rc c n ,rc c n rc c n 4^+^-][5-5+] JS ,h][S ,] 45_/4[5-5+] 

= /l[5 +5 +] - /i[5 -5 +] - Qr[5 -5 +] - (TsP -—-- (Ts cr,p -- p- 


[5+] 


[5+] 


[5-] 


^[5_/+][5-5+] [5_54([5-5+] + [5_/+]) [5_5+]([5_5+] + [5-/+]) 

- CTlP --h a - tt—-. - o- 


[5-] 


[5-] 


[5-] 


45+5+] , ,, , , .,rc c n . 4^+^-][5+5+] , 45+/+][5+5+] , , [5_5+]([5_5+] + [5-/+]) 

= 2a[iy +] — 2A[S +S +] — 2crsP -j44-7- — 2crscriP -7^;—;-h 2Qr- 


dt 


[5+] 


[5+] 


[54 


( 4 ) 






























































References 


References 

[1] A. L. Lloyd, R. M. May, How vinases spread among computers and people. Science 292 (5520) (2001) 1316-1317. 

doi:10.1209/0295-5075/82/38004 

[2] L. Danon, A. R Ford, T. House, C. Jewell, M. J. Keeling, G. O. Roberts, J. V. Ross, M. C. Vernon, Networks and the epidemiology of 
infectious disease. Interdisciplinary Perspectives on Infectious Diseases 2011 (2011) 284909. 

[3] R. Pastor-Satorras, C. Castellano, P. Van Mieghem, A. Vespignani, Epidemic processes in complex networks, arXiv: 1408.2701 (2014). 

[4] J. Sanz, C.-Y. Xia, S. Meloni, Y. Moreno, Dynamics of interacting diseases, Phys. Rev. X 4 (2014) 041005. 

[5] M. E. J. Newman, Threshold effects for two pathogens spreading on a network, Phys. Rev. Lett. 95 (2005) 108701. 

[6] Y.-Y. Chen, F. Chen, D. Gunnell, P. S. F. Yip, The impact of media reporting on the emergence of charcoal burning suicide in taiwan, PloS 

ONE8(2013)e55000. 

[7] M. De Domenico, A. Lima, P. Mougel, M. Musolesi, The anatomy of a scientific rumor, Sci. Rep. 3 (2013) 2980. 

[8] A. Montanari, A. Saberi, The spread of innovations in social networks, Proc. Natl. Acad. Sci. U.S.A. 107 (2010) 20196-20201. 
doi:10,1073/pnas,1004098107 

[9] D. Lazer, R. Kennedy, G. King, A. Vespignani, The parable of Google flu: traps in big data analysis 343 (2014) 1203-1205. 

[10] H. A. Cameiro, E. Mylonakis, Google trends: a web-based tool for real-time surveillance of disease outbreaks, Clin. Infect. Dis. 49 (2009) 
1557-1564. 

[11] World Health Organization. Consensus document on the epidemiology of severe acute respiratory syndrome (SARS), World Health Organi¬ 
zation, Geneva, Switzerland. 

[12] Z.-X. Tai, T. Sun, Media dependencies in a changing media environment: The case of the 2003 SARS epidemic in China, New Media Soc. 9 
(2007) 987-1009. 

[13] S. Funk, E. Gilad, C. Watkins, V. A. A. Jansen, The spread of awareness and its impact on epidemic outbreaks, Proc. Natl. Acad. Sci. U.S.A. 
106 (2009) 6872-6877. 

[14] S. Funk, M. Salathe, V A. A. Jansen, Modelling the influence of human behaviour on the spread of infectious diseases: a review, J. R. Soc. 
Interface 7 (2010) 1247-1256. 

[15] C. Granell, S. Gomez, A. Arenas, Dynamical interplay between awareness and epidemic spreading in multiplex networks, Phys. Rev. Lett. 
Ill (2013) 128701. 

[16] W. Wang, M. Tang, H. Yang, Y. Do, Y.-C. Lai, G. Lee, Asymmetrically interacting spreading dynamics on complex layered networks, Sci. 
Rep. 4 (2014)5097. 

[17] C. Granell, S. Gomez, A. Arenas, Competing spreading processes on multiplex networks: Awareness and epidemics, Phys. Rev. E 90 (2014) 
012808. 

[18] S. Funk, E. Gilad, V. A. A. Jansen, Endemic disease, awareness, and local behavioural response, J. Theor. Bio. 264 (2010) 501-509. 

[19] F. D. Sahneh, F. N. Chowdhury, C. M. Scoglio, On the existence of a threshold for preventive behavioral responses to suppress epidemic 
spreading, Sci. Rep. 2 (2012) 632. 

[20] H.-F. Zhang, Z.-X. Wu, M. Tang, Y.-C. Lai, Effects of behavioral response and vaccination policy on epidemic spreading - an approach based 
on evolutionaiy-game dynamics, Sci. Rep. 4 (2014) 5666. 

[21] Q.-C. Wu, X.-C. Fu, M. Small, X.-J. Xu, The impact of awareness on epidemic spreading in networks. Chaos 22 (2012) 013101. 

[22] A.-L. Barabasi, R. Albert, Emergence of scaling in random networks. Science 286 (1999) 509-512. 

doi:10,1126/science,286,5439,509 

[23] D. J. Watts, S. H. Strogatz, Collective dynamics in ‘small-world’ networks. Nature 393 (1998) 440^42. 

[24] P. Holme, J. Saramaki, Temporal networks, Phys. Rep. 519 (2012) 97-125. 

[25] R. Pastor-Satorras, A. Vespignani, Epidemic spreading in scale-free networks, Phys. Rev. Lett. 86 (2001) 3200. 

dci:10,1103/PhysRevLett,86,3200 

[26] C. Liu, Z.-K. Zhang, Information spreading on dynamic social networks, Commun Nonlinear Sci Numer Simula! 19 (4) (2014) 896-904. 

[27] C. Gracia-Lazaro, A. Ferrer, G. Ruiz, A. Taranc’on, J. A. Cuesta, A. S’anchez, Y. Moreno, Heterogeneous networks do not promote cooper¬ 
ation when humans play a Prisoner’s Dilemma, Proc. Natl. Acad. Sci. U.S.A. 109 (2012) 12922-12926. 

[28] T. Gross, C. J. D. D’Lima, B. Blasius, Epidemic dynamics on an adaptive network, Phys. Rev. Lett. 96 (2006) 208701. 

dci;10,1103/PhysRevLett,96,208701 

[29] T. Gross, B. Blasius, Adaptive coevolutionaiy networks: a review, J. R. Soc. Interface 5 (2008) 259-271. 

[30] T. Gross, H. Sayama, Adaptive networks, Springer, Berlin, 2009. 

[31] I. Tunc, L. B. Shaw, Effects of community structure on epidemic spread in an adaptive network, Phys. Rev. E 90 (2014) 022801. 

[32] Y-Z. Zhou, Y.-J. Xia, Epidemic spreading on weighted adaptive networks, Physica A 399 (2014) 16-23. 

[33] E. Volz, L. A. Meyers, Susceptible-infected-recovered epidemics in dynamic contact networks, Proc. R. Soc. B 274 (2007) 2925-2933. 

dci;10,1098/rspb,2007,1159 

[34] S. Wieland, T. Aquino, A. Nunes, The structure of coevolving infection networks, EPL (Europhys. Lett.) 97 (2012) 18003. 
dci:10,1209/0295-5075/97/18003 

[35] D.-C. Guo, S. Trajanovski, R. van de Bovenkamp, H.-J. Wang, P. Van Mieghem, Epidemic threshold and topological strticture of susceptible- 
infectious-susceptible epidemics in adaptive networks, Phys. Rev. E 88 (2013) 042802. 

[36] I. Tunc, M. S. Shkarayev, L. B. Shaw, Epidemics in adaptive social networks with temporary link deactivation, J. Stat. Phys. 151 (2013) 
355-366. 

[37] L. D. Valdez, P. A. Macri, L. A. Braunstein, Intermittent social distancing strategy for epidemic control, Phys. Rev. E 85 (2012) 036108. 

[38] A. J. Morris, Representing spatial interactions in simple ecological models, Ph.D. thesis. University of Warwick (1997). 

[39] M. J. Keeling, The effects of local spatial structure on epidemiological invasions, Proc. R. Soc. London, Ser. B 266 (1999) 859-867. 


8 


[40] N. Crokidakis, S. M. D. Queiros, Probing into the effectiveness of self-isolation policies in epidemic control, J. Stat. Mech. (2012) P06003. 

[41] C. Lagorio, M. Dickison, F. Vazquez, L. A. Braunstein, P. A. Macri, M. V. Migueles, S. Havlin, H. E. Stanley, Quai'antine-generated phase 
transition in epidemic spreading, Phys. Rev. E 83 (2011) 026102. 

[42] L. Lti, D.-B. Chen, T. Zhou, The small world yields the most effective information spreading. New J. Phys. 13 (2011) 123005. 
doi:10,1088/1367-2630/13/12/123005 

[43] Z.-K. Zhang, C.-X. Zhang, X.-P. Han, C. Liu, Emergence of blind areas in information spreading, PLoS ONE 9 (4) (2014) e95785. 


9 



Figure 1: (Color online) Transmission diagram of SI epidemic model and the information about the epidemic among the population. 
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Figure 2: (Color online) The epidemic spreading dynamics of vaiious information diffusion probabilities a without considering the effect of 
adaptive process. The parameters are set as = 0.3, crs = 0.5, cr/ = 0.7, A = 0.2, u = 0.75, rs = rj = 0. The inset shows the information diffusion 
dynamics {Info) of various for a = 0.6. 
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Figure 3: (Color online) Compaiison of simulation results with the mean-field model and the pair approximation model without considering the 
effect of adaptive process. The parameters are set as (5 = 0.3, crs = 0.5, (Ji - 0.7, /i = 0.2, o} = 0.75, a = 0.6, rs = rj = 0. 




T k 

Figure 4: (Color online) Dynamical analysis of the spreading model with adaptive process, (a) Comparison of the pairwise model with the 
simulation results, (b) Degree distribution of the original network and that after the adaptive process. The pai'ameters ai‘e set as = 0.2, cr^ = 
0.5,cr/ = 0.7, A = 0.2,0) = 0.2, a = 0.5, rs =0.\5,rj = 0.1. 
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Figure 5: (Color online) The fraction of infected individuals in the stationary state (colors in the phase diagram represent the density of infected 
individuals at the final state, the dashed green curve shows that the prevalence value transmits from near 0 to significantly larger than 0) versus a 
and yS for (a) pairwise analysis and (b) simulation result. The parameters are set as crs = 0.5, ctj = 0.7, A = 0.2, lj = 0.2, rs = r/ = r = 0.1. 



Figure 6: (Color online) Illustration to dynamic spreading process by considering the adaptive effect on the lattice. The square gridding patterns 
show the distribution of the infected and informed individuals in some particulai* time steps. The red ai'ea represents the nodes that are infected by 
the epidemic, while the gray area represents the informed individuals. The red curves (lower panels) describe the fraction of infected individuals 
over time with coiTesponding. (a) cr = 0.1; (b) a = 0.3; (c) a = 0.5; (d) a = 0.7. Other pai'ameters are set as yS = 0.4, (Ts = 0.4, (Tj = 0.8, d = 
0.1, tu = 0.2, rs = r] = r = 0.1. 
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Figure 7: (Color online) Fraction of infected individuals versus r. Different curves correspond to different a. Other parameters are set as yS = 
0.2, crs = 0.5, cr, = 0.7, A = 0.2, co = 0.2. 
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Figure 8: (Color online) Distribution of the infected density in the final state versus different values of r and /?. Each distribution is obtained by 
carrying out 10,000 independent realizations for the final fraction of infected. The parameters are set as = n = r = 0.7,0.35,0.15,0;yS = 
0.05,0.35,0.25,0.4 for (a), (b), (c) and (d) respectively. Other parameters are (Ts = 0.5, ct/ = 0.7, A = 0.2, cj = 0.2, a = 0.6. 



Figure 9: (Color online) (a) Bifurcation diagram of the density of the infected I as a function of the infection probability jA for different values of 
the edge-breaking rate r based on the results of simulation of the full network (diamonds), (b) Two parameter bifurcation diagram showing the 
dependence on the edge-breaking rate r and the infection probability (A based on the results of simulation of the full network, (c) Full phase diagram 
r — jA for the simulation of the adaptive process. The pai'ameters are the same as Fig. 8. 
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